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We investigated the Abrikosov vortex lattice (VL) of a pure Niobium single crystal with the muon 
spin rotation (/xSR) technique. Analysis of the pSR data in the framework of the BCS-Gor'kov theory 
allowed us to determine microscopic parameters and the limitations of the theory. With decreasing 
temperature the field variation around the vortex cores deviates substantially from the predictions 
of the Ginzburg-Landau theory and adopts a pronounced conical shape. This is evidence of partial 
diffraction of Cooper pairs on the VL predicted by Delrieu for clean superconductors. 
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The Ginzburg-Landau (GL) theory for superconduc- 
tors is expressed in terms of an order parameter A(r). 
While the absolute value of A(r) determines the local 
superfluid density, its phase gradient is proportional to 
the local magnetic vector potential. This phase variation 
leads to the magnetic flux quantization and the formation 
of a periodic vortex lattice (VL) in type-II superconduc- 
tors as was predicted by Abrikosov [l| . The VL field vari- 
ation is uniquely characterized by two length scales, the 
magnetic penetration depth A and the coherence length 
^GL • This simple model turned out to be quite successful 
in describing the behavior of a superconductor in a mag- 
netic field [2, 0] and serves as a basis for data analysis 
of experiments 0, Q. However, as was shown theoreti- 
cally by Delrieu [6|, the GL model is unable to describe 
the magnetic response in clean superconductors at low 
temperatures and close to the upper critical field i?c2- 

Soon after the publication of the microscopic Bardeen- 
Cooper-SchriefFer (BCS) theory for conventional super- 
conductors j7jj using a Green function formalism Gor'kov 
derived the GL equations from the BCS theory [8] . Based 
on Gor'kov's equation Delrieu analyzed the field varia- 
tion for classical s— wave superconductors in the vicinity 
of Bc2 i9]- He found that for clean superconductors the 
Cooper pairs (CPs) with balistic trajectories through the 
vortex cores diffract on the periodic potential of the VL. 
As a result, in the low temperature limit close to Bc2 the 
spatial field variation around a vortex core has a coni- 
cal shape, rather than the cosine-like GL behaviour, and 
the fields at the minimum and saddle points are interex- 
changed relative to the GL prediction. Nearly at the 
same time Brandt came to the same conclusion based 
on a nonlocal theory of superconductivity Q. To ob- 
serve the effect of diffraction of CPs, the carrier mean- 
free path ^mfp should exceed the intervortex distance, the 
measurements should be done in the vicinity of Bc2iT) 
such that the A(r) gradient is negligible, and the tem- 



perature should be low to minimize thermal fluctuations. 

Although the theoretical study of the influence of the 
diffraction of CPs on the field variation was already per- 
formed in 1972, it was not investigated experimentally 
in detail so far. Early muon spin rotation (/iSR) experi- 
ments revealed a linear high- field tail in the magnetic field 
distribution D'^^{B^) [10|, in agreement with the the- 
oretical expectation. However, as noticed recently [ll|, 
it occurred at an unexpectedly high temperature. As we 
note below the temperature stability is crucial in order to 
minimize experimental artifacts also leading to a linear 
high-field tail in D'^''p{B^). On the other hand, a /xSR 
study of vanadium did not reveal any deviation from the 
GL theory [l^. Thus, a superconductivity in the clean 
limit is one of the critical conditions for the observation of 
the high-field linear tail. Most of the novel high temper- 
ature superconductors (HTSs) are in the clean limit, and 
the tail should be observed provided the measurements 
are performed close to Bc2 and at low temperature. Such 
studies of HTSs still await to be performed. 

As a first superconductor to look for the effect of 
CP diffraction, we have chosen metallic niobium (Nb), 
since it is a simple BCS superconductor and pure sin- 
gle crystals are available. It is a type H superconductor 
(k « 0.8 > ~ 0.7), and therefore a VL is expected 

in the bulk when an external field Boxt larger than the 
lower critical field is applied. As shown by small an- 
gle neutron scattering, for Bext parallel to the crystallo- 
graphic (111) direction the VL exhibits a simple triangu- 
lar lattice 13- 

Our Nb sample was a single crystal disk of 13 mm 
diameter and 2 mm thickness with a (111) axis oriented 
normal to the disk. The samples studied here and in 
Ref. [l^l come from the same batch, so they should be 
of the same metallurgical quality. The sample is further 
characterized when discussing its value of -Bc2(0). 

The /xSR experiments were performed at the Swiss 
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Muon Source (S/LtS), Paul Scherrer Institute (PSI), 
Switzerland, using the general purpose spectrometer 
(GPS) for T > 1.6 K and the low temperature facility 
(LTF) for T < 1.6 K. A field cooled procedure was used 
with Bext perpendicular to the sample plane (parallel to 
a (111) axis). The ^SR spectra were recorded with the 
transverse field geometry, i.e. the initial muon spin polar- 
ization was perpendicular to Bcxt- By definition, Bext 
is parallel to the Z axis of the laboratory orthogonal ref- 
erence frame. With this geometry the field distribution 
in the bulk of a superconductor can be probed [3l ■ Since 
Bcxt was much larger than the field distribution width, it 
is its component along the Z axis which was measured, 
i.e. Dl''^{B^). We explicitly distinguish Dc{B^) from 
De^P(i3^), since Dc{B^) only stands for the distribution 
of a perfect VL without crystal disorder. 

The forward and backward positron detectors with re- 
spect to were used to build the /zSR asymmetry time 
spectra A{t) recorded with total statistics ranging from 
1.0 X 10'' to 8.0 X 10^ positron events. Typical A{t) in 
the normal and the superconducting states are displayed 
in the insert (a) of Fig. [T] Note that in contrast to the 
normal state, a strong damping of A{t) in the supercon- 
ducting state is observed which is characteristic of the 
local magnetic field variation due to the VL. From these 
kind of measurements Bc2{T) was determined, yielding 
Be2(0) = 430(2) mT (see Fig. [T]). This value is smaller 
than -Bc2(0) = 443 mT reported for a sample with a 
residual resistivity ratio RRR = 750 [l^l ■ Hence, for our 
sample the RRR > 750. Our value of i3c2(0) indicates 
that the sample is of high quality and pure [3]. The 
analysis of the /iSR spectra below supports this conclu- 
sion. 

As shown in the insert (b) of Fig. [1] the temperature 
stability is important for recording high quality data close 
to Bc2- Large fluctuations of temperature may lead to 
substantial smearing of the measured spectra. The ex- 
perimental and theoretical field distributions presented 
in this letter, were obtained by Fourier transformation 
(FT) of the Gaussian apodized time spectra (i.e. Fourier 
transform of A(t) exp(— cr^ppt^/2) where Capp = 4.7 /zs). 
Note that the apodization has no influence on the anal- 
ysis, smce we directly fit A{t), rather than D^''p(B^). 

Figure [2] displays the field distributions measured in 
the vicinity of Bc2{T) from a temperature of 7.8 K close 
to Tco (critical temperature at low field) down to 0.02 K 
[see Fig. [T] for the location of the points (crosses) in the 
field/ temperature diagram]. For each D'^'^{B^) mea- 
sured at LTF a relatively intense sharp peak is present 
at a field slightly larger than Bcxt, in contrast to the 
GPS data for which only a small hump is found. This 
field structure (sharp peak and small hump) originates 
from the muons stopped in the cryostat walls and sam- 
ple holder (background signal). The present GPS data 
are therefore of much better quality than previous results 
[lo| obtained in the same temperature range showing an 
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FIG. 1. (Color online) Bc2(r) for Bcxt || (111) as deter- 
mined by /iSR measurements for our Nb single crystal sam- 
ple (circles). The dashed line corresponds to the equation 
Sc2(r) = Bc2(0)(l -r2)/(l + r2) proposed in Ref. [TJ. Here 
r = r/Tco with Tco = 9.25 K and Bc2(0) = 430 (2) mT. The 
crosses indicate the points in the field-temperature diagram 
at which the field distributions displayed in Fig. [2] were mea- 
sured. Insert (a) shows /^SR asymmetry spectra A(t) recorded 
at 1.6 K in the normal (o) and in the mixed (•) states for 
Bcxt = 450 and 360 mT, respectively. The solid lines are fits 
of Eq. ([l]) to the data. The spectra are shown in a rotat- 
ing frame of 440 and 350 mT, respectively. Insert (b) shows 
the effect of the sample temperature stability on D'^^'^(B^). 
The two measurements were performed at _Bcxt = 165 mT and 
T — 5.9 K with different temperature regulation systems. For 
one of them the temperature was found to oscillate periodi- 
cally around an average with a period of few seconds and a 
peak to peak amplitude of 150 mK. For the other system this 
amplitude was reduced to about 20 mK. 

intense background signal. Qualitatively, a linear high- 
field tail in Dc^p{B^) is inferred at maybe 1.2 K and 
certainly at the lower temperatures, but not at higher 
temperature. On the other hand, this tail is already seen 
at 2.6 K in the pubhshed data 

A /iSR spectrum is described by the sum of two con- 
tributions: 

A{t) = Ao ■ [F,Rs{t) + (1 - F,)i?bg(t)], (1) 

where Aq is the total initial asymmetry and Fs the frac- 
tion of muons stopped in the sample. We have Aq — 
0.211 (2) [0.217 (2)] and = 0.993 (2) [0.817 (2)] for mea- 
surements carried out with the GPS [LTF] spectrometer. 
The function 

Rs{t) = e-*"'*' J Dc{B^) cos(7^B^t + M^B^ , (2) 
describes the time evolution of in the sample while 
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FIG. 2. (Color online) Field distributions D^'^'^iB^) of Nb single crystal obtained in the vicinity of Bc2(T) with the LTF 
spectrometer (a) and the GPS (b). The solid lines represent best fits of Eqs. ((Ij-Q to the data. As explained in the text, 
6 = 0.01 and a is a free parameter for all the asymmetries A{t). Leaving the last available parameter vp free, leads to proper 
fits for T > 0.8 K, but not below. The blue dashed-dotted lines in panel (a) represent the results for T — 0.5, 0.1, and 0.02 K. 
At 0.5 K the misfits are small, but still present. The solid lines for T < 0.5 K are computed with the fixed value c = 0.08 
for T = 0.6 K. The blue dotted lines in panel (b) correspond to fits with NGL, while the red dashed and blue dashed-dotted 
lines visualize the differences between Df^'^{B^) and the BCS-Gor'kov and NGL predictions, respectively. Panels (c) and (d) 
illustrate field variation and contour plot of (v) at 7.8 and 0.5 K, with denoting the field at the vortex core. 



Rhgit) — exp(— cos(7^i3^gt + (/)o) accounts for 
the background. Here 7^ = 851.6 Mrads^^T^^ is the 
muon gyromagnetic ratio, ahg — 0.22(2) /is~^ stands for 
the background damping, and 0o is the initial phase. The 
mean field for the background B^^ is only slightly differ- 
ent from Bcxt (see Fig. [2]). We write cr'^ = a^^ + cr^jg, 
where a^u accounts for the damping due to the nuclear 
^■^Nb spins, and ddis is the parameter describing the ef- 
fect of VL disorder. As usual, we assume the effect of 
disorder to be modelled by a Gaussian function with a 
field standard deviation Udia/l^i. ^2Q] . Although this is a 
crude approximation, we note that the influence of VL 
disorder on the high- field tail of D'^'^{B^) is relati vely 
moderate compared to its effect on the low- field side [21 1 . 



We determine Dc{B^) from the real space field 
map B^ [y) of the two-dimensional VL: Dc{B^) = 
J^^ 6{B^{r) — B^)(Pr, where the integral extends over 



the VL unit cell. In terms of its Fourier components, 

B^ir) = ^l,„,. exp(zK,„,, • r), (3) 

where the sum is over the reciprocal space. 

First we analyze the data with numerical solution of 
the GL (NGL) model for B^ ^ using Brandt's method 
01 . It only depends on A and ^gl- The fit for the high- 
est temperature data is reasonable (see Fig. [^b)- We 
get A = 59.4 (2) nm and as = 0.72 (1) mT, with ^gl = 
66.5 (2) nm estimated from the measured i3c2(7.8 K) and 
using the GL formula: ^gl = ($o/27ri?c2)^^^- Here, $0 is 
the flux quantum. As expected, k = A/^gl = 0.89 (1) > 
I/V2 ~ 0.7. The results for T < 5.9 K in Fig. ^ were 
obtained with A, ^gl, and as as free parameters. The GL 
model fails to describe the high- field tails in D°^p{B^). 
In addition, unreasonable large k values are derived. For 
example: k = 48.7/28.5 — 1.7 at 1.6 K. If k had been 
taken temperature independent as it should, the misfits 
are even worse. As expected, the GL model can only 
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describe D^'^p(B^) well near Tco. 

Next we analyze the data with the BCS-Gor'kov the- 
11, [23] ■ First we discuss the characteristics of 
in the vicinity of Bc2{T). We use the notations 
The Fourier component is a function of four 
parameters: B§_^ ^ — /^^^(a, &, c, J) [1, [HI, [ll]- Here, 

a — —iiottNoAq2c/ does not influence the shape of 
B^ (r) and therefore Dc{B^), but only determines the 
scale of the field variation. It is proportional to the den- 
sity of states at the Fermi level Nq (per spin, volume, 
and energy), the quantity Aq = |A(r)p (|A(r)P is the 
spacial average of |A(r)p), and is inversely proportional 
to the spatial average of the field B^{r). The dimen- 
sionless parameters b, c, and d determine the shape of 
Dc{B^) and are expressed by the ratios of four length 
scales: b = {A/tt£_^)^, c = A/£_'^, and d = A/£nifp- 
Here, A = [$o/(27ri3-^)]^/^ is a length parameter propor- 
tional to the intervortex distance. The field and temper- 
ature dependent length scale — fivp /{Tf^o) diverges 
a.tW -i' Sc2, i.e. b 0, while = hvF/{27rkBT). 
The parameter c is strongly temperature and field depen- 
dent. It vanishes as T — and diverges at T — >• Tco [11| . 
Finally, for clean superconductors d is negligibly small, 
since imfp significantly exceeds the intervortex distance. 

Cooper pair diffraction may influence Dc{B^) when 
three experimental conditions are met: A <^ imfp, A <^ 
ttC^, and A < HjIIII,!!!. The first condition implies 
a clean superconductor, the second is only satisfied in the 
vicinity of Bc2{T), and the third one is only possible at 
low T. Thus, the minimum of {imip, t^^^ , } deter- 
mines the effective diffraction length scale of CPs relative 
to the intervortex distance 2.693 x A. 

Assuming simple conventional temperature depen- 
dences for Aq and ^gl, and in addition a standard field 
dependence for Ag, 6 is computed after selecting values 
for Tco and -Bc2(0) 0. For Nb we get 0.01 < & < 0.02, 
except for the spectrum at 7.8 K for which b = 0.04. Re- 
calling that B^ is only weakly b dependent for such 
small values, the approximate nature of the field and 
temperature dependences we have taken for b is not of 
importance. The data analysis was done with b = 0.01. 
Since 
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the Fourier component B^ ^ depends on a, d and vp. 

The analysis of A{t) for T > 0.8 K shows that d < 
0.01, which agrees with the estimate of fmfp — 7 /xm for 
the Nb sample with RRR= 750 [11]. Consequently, we 
are in the clean limit and d has a negligible influence 
on Dc{B^). The results of the analysis are presented 
in Fig. [2j The BCS-Gor'kov model describes the high- 
field tail significantly better at low temperatures while at 
the highest temperature both models reproduce the data 



equally well. The deviation from the GL theory and the 
gradual disappearance of the cut-off singularity at the 
maximal field is a result of the conical shape of the spatial 
field variation at the vortex cores (see Fig.[2ji), which in 
turn is a consequence of the partial CP pair diffraction. 
The BCS-Gor'kov model breaks down for T < 0.5 K as 
shown by the dashed-dotted lines in Fig. [2^. A proper 
description requires to consider B^ as a function of 
d and c rather than of a and vp, and to keep the value 
of c at T = 0.6 K for the lower temperatures (solid lines 
in Fig. [5^). This means that the sharpness of the B^{r) 
cones is limited by a physical process. Refering to Eq. ([3]), 
we suggest that the VL structural disorder may round up 
the cones, as observed experimentally. 

We get VF = 2.0(2) x 10^ m/s from the fits for 
T > 0.8 K. This value is smaller but compatible with 
i;f 2.73 X 10^^ [H and 2.94 x 10^ m/s [H determined 
from magnetization measurements. From the measured 
d and c we can compute —B^d/c ~ 2ttiiqNq/\q. While 
iVo is a constant, Aq is field and temperature dependent. 
At T = and interpolating Aq to zero field with a con- 
ventional formula [Ilj, [-Wd/c]/[l - [B^ / B^2{Q))] = 
2t:^oNo/!\q{Q) where Ao(0) is taken as the BCS gap for 
s— wave superconductors which is temperature indepen- 
dent for T ^- [2j. From Ao(0) = 1.45 meV [I^ and 
our estimate 27r/xoiVoAg(0) = 0.0155(5) T^, we obtain 
No w 3.4 X 10^6 J-im-3spin-i w 0.20 eV-V-ispin-\ 
where v^.c. stands for the volume of the Nb body cen- 
tered cubic unit cell. This can be compared to the spe- 
cific heat result A^o = 0.42 eV~^v-i spin'^ [24]. Al- 
though smaller by a factor 2, our value is still reason- 
able recalling the approximate nature of the interpolat- 
ing field formula. In our measurements all the conditions 
for the observation of partial CP diffraction are met at 
T < 1.2 K: A/tt^^ < 0.11 < 1, A/Cfp < 0.01 < 1, and 
A/^'^ = 0.16 < 1. 

To conclude, we investigated the magnetic field distri- 
bution for the vortex lattice of a pure Nb single crystal 
with the ^SR technique. The data were analyzed us- 
ing the solution of the BCS-Gor'kov equation proposed 
by Delrieu [^. As a result, we found strong evidence 
for partial Cooper pair (CP) diffraction on the periodic 
potential of the vortex lattice reflected in the conical nar- 
rowing of the real space field variation around the vortex 
cores. However, the zero-temperature limit prediction is 
not realized experimentally, presumably due to the resid- 
ual disorder of the vortex lattice. From the analysis we 
determined two microscopic parameters: the Fermi ve- 
locity Up and the product of the density of states at the 
Fermi level A^o and square of the order parameter at zero 
temperature Aq(0). The signature of partial CP diffrac- 
tion in small angle neutron scattering data is the presence 
of many high order Bragg's reflections in the vicinity of 
Bc2{T) with the predicted intensity ratios 
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This is 

an experimental challenge worthwhile to attempt. The 
observation of partial CP diffraction should not be re- 
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stricted to Nb. Under proper experimental conditions it 
should also be seen for any clean type-II superconductor. 

This work was performed at the Swiss Muon Source 
(S/iS), Paul Scherrer Institut (PSI), Switzerland and 
partly supported by NCCR MaNEP sponsored by the 
Swiss National Science Foundation. 



1174 (1957). 

Superconductivity 



[1] A. A. Abrikosov, Sov. Phys. JETP 5, 
[2] M. Tinkham, Introduction to 

(McGraw-Hill, New York, 1996). 
[3] J. B. Ketterson and S. N. Song, Superconductivity (Cam- 
bridge University Press, Cambridge, 1999). 
E. H. Brandt, Phys. Rev. Lett. 78, 2208 (1997), URL 
http : //link. aps . org/doi/10 . 1103/PhysRevLett ^78 . 2208 
J. E. Sonier, Reports Prog. in 

Phys. 70^ 1717 (2007), URL 

|http : //iopscience . iop . org/0034-4885/70/ll/R01 | 

[6] J. M. Delrieu, J. Low Temp. Phys. 6, 197 (1972), URL 

http:/7wwwTspringerlink."cWcontent/xvl26il085wj3vp6/?MUB^Fliel, Fiz. Met. Metalloved. 37, 63 (1974) 



[4] 
[5] 



J. E. S onier, Phys. Rev. B 74, 054511 (2006), URL 
lhttp://link . aps . org/doi/10 . 1103/PhysRevB . 74 . 0545lT] 
[13] J. Schelten, H. UUmaier, and W. Schmatz, physica sta- 
tus solidi (b) 48, 619 (1971), ISSN 1521-3951, URL 
http : //dx . doi . org/10 . 1002/pssb . 2220480219 
[14] R. Kahn and G. Parette, Solid State Com- 

mun. 13, 1839 (1973), ISSN 0038-1098, JJRL 

http: //www. sciencedirect . com/science/cirticle/pii/003810i 
[15] S. Miihlbauer^ C. Pfleiderer, P. Boni, M. Laver, 
E. M. Forgan, D. Fort, U. Keiderling, and 

G. Behr, Phys. Rev. Lett. 102, 136408 (2009), URL 

http : //link . aps . org/doi/10 . 1 103/PhysRevLett . 102. 13640'8] 
[16] A. Yaouanc and P. Dalmas de Reotier, Muon Spin Rota- 
tion, Relaxation, and Resonance: Applications to Con- 
densed Matter, International Series of Monographs on 
Physiscs 147 (Oxford University Press, Oxford, 2011). 
[17] D. K. Finnemore, T. F. Stromberg, and C. A. 
Swenson, Phys. Rev. 149, 231 (1966), URL 
http://link.aps.org/doi/10. 1103/PhysRev. 149.231 
[18] S. J. W illiamson, Phys. Rev. B 2, 3545 (1970), URL 

|http:/ /link . aps . or g/doi/10 . 1103/PhysRevB . 2 . 3545 
[19] N. Y. Alekseyevskiy, V. I. Nizhankovskiy, and H.-H. 



[7] J. Bardeen, L. N. Cooper, and J. R. Schri- 
effer, Phys. Rev. 108, 1175 (1957), URL 
http://link.aps.org/doi/10. 1103/PhysRev. 108. 1175 
[8] L. Gor'kov, Sov. Phys. JETP 36, 1364 (1959)"^ 
[9] E. H. Brandt, Phys. Stat. Sol. B 
64, 467 (1974), ISSN 1521-3951, URL 
http :7/dx . doi . org/10 . 1002/pssb . 22210640207 
[10] D. Herlach, G. Majer, J. Major, J. Rosenkranz, 
M. Schmolz, W. Schwarz, A. Seeger, 
W. Tempi, E. Brandt, U.Essmann, et al., 
Hyperfine Interactions 63, 41 (1990), URL 



[20] E. H. Brandt, J. Low Temp. Phys. 73, 355 (1988), URL 

http://www.springerlink.com/content/p8442442n7486186/7ffl 

[21] A. Maisuradze, R. Khasanov, A. Shen- 
gelaya, and H. Keller, J. Phys.: Con- 
dens. Matter 21, 075701 (2009), URL 
http : //stacks . iop . org/0953-8984/21/i=7/a=075701 

[22] J. M. Delrieu, Ph.D. thesis. University of Paris-Orsay 
(1974). 

[23] H. R. Kerchner, D. K. Christen, and S. T. 
Sekula, Phys. Rev. B 24, 1200 (1981), URL 



[11] 
[12] 



^tp : //link . aps . org/doi/10 . 1103/PhysRevB . 24 . 1200l 

http : //www . springerlink . com/ content /?^=puh7.3a.Va8nypefi4\n^±mmmt,%Q^Vimm^ 



A. Maisuradze and A. Yaouanc, arXiv:1302.5041 (2013), 

URL http : //arxiv . org/abs/1302 .5041 

M. Laulajainen, F. D. Callaghan, C. V. Kaiser, and 



Temperature Physics 62, 301 (1986), ISSN 0022-2291, 
URL http : //dx ■ doi . org/10 ■ 1007/BF00683466 1 



